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Beschreibung 

[0001] Die Erfindung betrifft ein Verfahren zum modellgestutzten Bestimmen der in die Zylinder einer Brennkraftma- 
schine einstromenden Luftmasse nach dem Oberbegriff des Patentanspruches 1. 

s [0002] Motorsteuerungssystemetur Brennkraftmaschinen, die mit Kraftstoffeinspritzung arbeiten, benotigen die vom 
Motor angesaugte Luftmasse m^^^als ein MaB fur die Motorlast, Diese GroBe bildet die Basis zur Realisierung eines 
geforderten Kraftstoff-Luft-Verhaltnisses. Wachsende Anforderungen an Motorsteuerungssysteme, wie die \ferrlnge- 
rung der Schadstoffemission von Kraftfahrzeugen, bedingen, daB die LastgroBe fur stationare und instationare Vor- 
gange mit geringen zuiassigen Fehlern bestimmt werden muB. Neben den genannten Betriebsfallen bietet die genaue 

10 Lasterfassung wahrend der Warmlauf phase der Brennkraftmaschine ein erhebliches Potential zur Schadstoffreduktion. 
[0003] Bei luftmassengefuhrten Motorsteuerungssystemen stellt im Instationarbetrieb das als l^stsignal der Brenn- 
kraftmaschine dienende Signal des Luftmassenmessers, der stromaufwarts des Saugrohrs angeordnet ist, kein MaB 
fur die tatsachliche Fullung der Zylinder dar, weil das Volumen des Saugrohrs stromabwarts der Drosselklappe als 
Luftspeicher wirkt, der bef Dllt und entleert werden muB. Die maBgebende Luftmasse fur die Einspritzzeitberechnung 

IS ist aber diejenige Luftmasse, die aus dem Saugrohr heraus und in den jeweiligen Zylinder hineinstromt. 

[0004] Bei saugrohrdruckgeluhrten Motorsteuerungssystemen gibt zwar das Ausgangssignal des Drucksensors die 
tatsachlichen Druckverhaltnisse im Saugrohr wieder, die MeBgrbBen stehen aber u.a. auf grund der notwendlgen Mit- 
telung der MeBgroBe erst relativ spat zur Verfugung. 

[0005] Mit der Einf uhrung variabler Ansaugsysteme und variabler Ventilsteuerungen entstehen fur empirisch gewon- 
20 nene Modelle zur Gewinnung der LastgroBe aus MeBsignalen eine sehr groBe Vielzahl von EinfluBgroBen, die die 

entsprechenden Modellparameter beemflussen. 

[0006] Auf physikalischen Ansatzen basierende model Igestutzte Berechnungsmethoden stellen einen guten Aus- 
gangspunkt zur genauen Bestimmung der Luftmasse m2y/dar. 

[0007] Aus der DE 39 1 9 488 C2 ist eme Vorrichtung zur Regelung und zur Vorausbestimmung der Ansaugluftmenge 
2S einer saugrohrdruckgefuhrtcn Brennkraftmaschine bekannt, bei der der Drosselklappenoffnungsgrad und die Motor- 
drehzahl als Grundlage zur Berechnung des derzeitigen Wertes der in den Brennraum der Maschine eingesaugton 
Luft verwendet werden. Diese berechnele. gegenwartige Ansaugluftmenge wird dann als Grundlage zur Berechnung 
des vorausbestimmten Wertes fur die Ansaugluftmenge, die in den Brennraum der Maschine zu einer bestimmten Zeit 
von dem Punkt an. an dem die Berechnung ausgefuhrt wurde, einzusaugen Ist, benutzt. Das Drucksignal, das strom- 
30 abwarts der Drosselklappe gemessen wird, wird mit Hilfe von theoretischen Beziehungen korrigiert, so daB eine Ver- 
besserung der Bestimmung der angesaugten Luftmasse erreicht und damit eine genauere Berechnung der Einspritz- 
zeit moglich isL 

[0008] Im instationaren Betrieb der Brennkraftmaschine ist es aber wunschenswert, die Bestimmung der in die Zy- 
linder einstromenden Luftmasse noch genauer durchzufuhren. 
35 [0009] Der Erfindung liegt die Aufgabe zugrunde ein Verfahren anzugeben, mit dem die tatsachlich in den Zylinder 
der Brennkraftmaschine einstromende Luftmasse mit hoher Genauigkelt bestimmt werden kann. AuBerdem sollen 
systembedingte Totzeiten. die aufgrund der Kraftstoffvorlagerung und der Rechenzeit bei der Berechnung der Ein- 
spritzzeit auftreten konnen, kompensiert werden. 

[0010] Diese Aufgabe wird gemaB den Merkmalen des Patentanspruches 1 gebst 
40 [0011] Vorteilhafte Weiterbildungen finden sich in den Unteranspruchen. 

[0012] Ausgehend von einem bekannten Ansatz ergibt sich eine Modellbeschreibung, die auf einer nichtlinearen 
Differentialgleichung basiert. Im folgenden wird eine Approximation dieser nichtlinearen Gleichung vorgestellL Im Er- 
gebnis dieser Approximation laBt sich das Systemverhalten mittels einer bilinearen Gleichung beschreiben, die die 
schnelle Losung der Beziehung im Motorsteuergerat des Krattfahrzeugs unter Echtzeitbedingungen gestattet. Der 
45 gewahlte Modellansatz beinhaltet dabei die Modellierung von variablen Saugsystemen und Systemen mit variablen 
Ventilsteuerungen. DIedurch diese Anordnung und durch dynamische Nachladung. d.h. durch Reflexionen von Druck- 
wellen im Saugrohr hervorgerufenen Effekte, konnen ausschlieBlich durch die Wahl stationar bestimmbarer Parameter 
des Modelles sehr gut berucksichligt werden. Alle Modellparameter sind einerseits physikalisch interpretierbar und 
andererseits ausschlieBlich aus stationaren Messungen zu gewinnen. 

[0013] Die meisten Algorithmen zur zeitdiskreten Losung der Differentialgleichung, die das Verhalten des hier ge- 
nutzten Modelles beschreibt. erfordern vor allem bei geringem Druckabfall uber der Drosselklappe, d.h. bei Vollast 
eine sehr kleine Rechenschrittweite, um numerisch stabil zu arbeiten. Die Folge ware ein unvertretbarer Rechenauf- 
wand bei der Bestimmung der LastgroBe. Da Lasterfassungssysteme meist segmentsynchron arbeiten, d.h. fur 4-Zy- 
lindormotoren wird alle 180° KW ein MeBwert abgetastet, muB die Modellgleichung ebenfalls segmentsynchron gelost 
werden. Im nachfolgenden wird ein absolut stabiles Differenzenschema zur Losung von Differentialgieichungen ein- 
gesetzt, das numerische Stabilitat bei belieblger Schrittweite garantiert. 

[0014] Das erfindungsgemaBe model Igestutzte Berechnungsverfahren bietet zudem die Moglichkeit einer Pradiktion 
des Lastsrgnales um eine wahlbare Anzahl von Abtastschritten, d.h. eine Vorhersage des Lastsignales mit variablem 
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Pradiktionshorizont. Wird die dem Pradiktionshorizont bei konstanter Drehzahl proportionale Pradiktionszelt nicht zu 

gro3, so erhalt man ein pradiziertes Lastsignal hoher Genauigkeit. 

[0015] Eine solche Vorhersage ist notwendig, da zwischen der Erfassung relevanter MeQwerte und der Berechnung 
der LastgroGe eine Totzeit entstelit. Deswelteren muB aus Grunden der Gemischaufbereitung vor dem eigentlichen 

s Beginn der Ansaugphase des jeweiligen Zylinders moglichst genau die Kraftsloffmasse uber die Einspritzventile zu- 
gemessen werden, die im Verlauf der kommenden Ansaugphase im gewunschten Verhaltnis zur Luftmasse m^^steht. 
Ein variabler Pradiktionshorizont verbessert die Gute der Kraftstoffzumessung im Instattonaren Motorbetrieb. Da bei 
steigender Drehzahl die Segmenlzeit abnimmt, mu3 der Einspritzvorgang eine groBere Anzahl von Segmenlen eher 
beglnnen, als dies bei einer niedrigeren Drehzahl der Fall ist. Um die zu doslerende Kraftstoffmasse moglichst exakt 

10 bestimmen zu konnen, ist die Pradiktion der Lastgr63e um die Anzahl von Segmenten^ um die die Kraftstoffvorlagerung 
vorgenommen wird, notwendig, um ein gefordertes Kraftstoff-Luft-Verhaltnis auch in diesem Fall einzuhalten. Die Pra- 
diktion der LastgroBe tragt somit aus einer wesentlichen Verbesserung der Einhaltung des gefordertene Kraftstoff- 
Luft-Verhaltnisses im instatlonaren Motorbetrieb bei. Dieses System zur modellgestutzten Lasterfassung ist in den 
bekannten Motorsteuerungssystemen, d.h. bei luftmassengefuhrte bzw. saugrohrdruckgefuhrte Motorsteuemngssy- 

15 steme wird im folgenden ein Korrekturalgorithmus in Form eines Modellregelkreises formuliert, der bei auftretenden 
Ungenauigkeiten von Model Iparametern eine permanente Genauigkeitsverbesserung, d.h. einen Modellabgleich im 
statk>naren und instatlonaren Betrieb gestattet. 

[0016] Ein Ausfuhrungsbeispiel des erfindungsgemaBen Vertahrens wird anhand der nachfolgenden schematischen 
Zeichnungen beschrieben. Dabei zeigen: 

20 

Figur 1 eine Prinzipskizze zum Saugsystem einer Otto-Brennkraftmaschine einschlieBlich der entsprechenden Mo- 
dell- und MeBgroBen, 
Figur 2 die DurchfluBfunktion und die dazugehorige Polygonzugapproximation, 

Figur 3 eine Prinzipdarstellung zum Modellregelkreis fur luftmassengefuhrte Motorsteuerungssysteme und 
2S Figur 4 eine Prinzipdarstellung zum Modellregelkreis fur saugrohrdruckgefuhrte Motorsteuerungssysteme. 

[0017] Bei der modellgestutzten Berechnung der LastgroBe m^y/ wird von der in Figur 1 dargestellten prinzipiellen 
Anordnung ausgegangen. Aus Grunden der Ubersichtlichkeit wird dabei nur ein Zylinder der Brennkraftmaschine dar- 
gestellt. Mit dem Bezugzeichen 10 ist dabei ein Saugrohr einer Brennkraftmaschine bezeichnet, in dem eine Drossel- 

30 klappe 11 angeordnet ist. Die Drosselkiappe 11 ist mit einem, den Offnungsgrad der Drosselklappe ermittelnden Dros- 
selklappenstellungsfuhler 14 verbunden. Stromauf warts der Drosselklappe 11 ist bei einem luftmassengefuhrten Mo- 
torsteuerungssystem ein Luftmassennnesser 1 2 angeordnet, wahrend bei einem saugrohrdruckgefuhrten Motorsteue- 
rungssystem ein Saugrohrdruckfuhler 1 3 im Saugrohr angeordnet ist. Je nach Art der Lasterfassung ist somit nur eine 
derbeidenKomponenten 12, 13vorhanden. Die Ausgange des Luftmassenmessers 12. des Drosselklappenstellungs- 

3S gebers 14 und des zum Luftmassenmesser 1 2 alternativ vortiandenen Saugrohrdrucksensors 1 3 sind mit EIngangen 
einer nicht dargestellten, ansich bekannten elektronischen Steuerungseinrichtung der Brennkraftmaschine verbunden. 
AuBerdem sind in Figur 1 noch schematisch ein EinlaBventil 15, ein AuslaBventil 16, sowie ein in einem Zylinder 17 
beweglichen Kolben 18 dargestellt. 

[0018] AuBerdem sind in Figur 1 ausgewahlte GroBen bzw. Parameter des Saugsystems eingezeichnet. Dabei be- 
40 deutet das Dachsymbol ""^^ uber einer GroBe. daB es sich um eine ModellgroBe handelt, wahrend GroBen ohne Dach- 
symbol MeBgroBen represent ieren. Im einzelnen bedeuten: 

[0019] P^ Umgebungsdruck, Pg Saugrohrdruck, Temperatur der Luft im Saugrohr. Vg das Volumen des Saug- 
rohrs. 

[0020] GroBen mit einem Punktsymbol kennzeichnen die erste zeitliche Ableitung der entsprechenden GroBen. m^)^ 
ist somit der Luftmassenstrom an der Drosselklappe und m^yi ist der Luftmassenstrom der tatsachlich In den Zylinder 
der Brennkraftmaschine einstromi. 

[0021] Die grundlegende Aufgabe bei der modellgestutzten Berechnung des Motorlastzustandes besteht nun in der 
Losung der Differentlalgleichung fur den Saugrohrdruck 



^DK " '"Zy/ 



(2.1) 



die sich unter der Voraussetzung konstanter Temperatur der Luft im Saugrohr T3 aus der Zustandsgleichung kdealer 
Gase herleiten laBt. 

[0022] Mit R|_ ist dabei die allgemeine Gaskonstante bezeichnet. 
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[0023] Die LastgroBe m^y/ wird durch Integration aus dem Zylindernnassenstronnl^v, bestimnnt. Die durch (2.1) be- 
schriebenen Verhaltnisse sind auf Mehrzylinder-Brennkraftmaschinen mit Schwingrohr- (Schaltsaugrohr-) und/oder 
Resonanzsaugsysteme ohne strukturelle Anderungen anwendbar. 

[0024] Fur Sysleme mit Multi-Point-Einspritzungen, bei denen die Kraftstoffzumessung durch mehrere Einspritzven- 
tile erfotgt. gibt die Gleichung (2.1) die Verhaltnisse genauer wieder als dies bei Single-Point-Einsprltzungen. d.h. bei 
Einspritzungen, bei denen der Kraftstoff mittels eines einzigen Kraftstoffeinspritzventiles zugemessen wird, der Fall 
ist. Bei erstgenannter Art der Kraftstoffzumessung 1st nahezu das gesamte Ansaugsystem mit Luft gefullt. Lediglich in 
einem kleinen Bereich vorden Einlaftventilen befindetsichein Kraftstoff-Luftgemisch. Im Gegensatz dazu ist bei Single- 
Point-Einspritzsystemen das gesamte Saugrohr von der Drosselklappe bis zum EinlaSventil mit Kraftstotf-Luft-Ge- 
misch gefullt, da das EInspritzventil vor der Drosselklappe angeordnet ist. In diesem Fall stellt die Annahme eines 
idealen Gases eine starkere Naherung dar, als dies bei der Multi-Point-Einspritzung der Fall ist. Bei Single-Point- 
Einspritzung erfolgt die Kraftstoffzumessung entsprechend^^^ bei Multi-Point-Einspritzung entsprechend^o^ 
[0025] Im folgenden wird die Berechnung der f^assenstrome OT£>A Und mzyi naher beschrieben. 
[0026] Die ModellgroBe des Luftmassenstromes an der Drosselklappe I^k wird durch die DurchfluBgleichung idealer 
Gase durch Drosselstellen beschrieben. An der Drosselstelle auftretende Stromungsverluste werden durch den redu- 
zierten Stromungsquerschnrtt A^^d berucksichtigt. Der Luftmassenstrom wird demnach durch die Beziehung 
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fur uberkritische Druckverhaltnisse bzw. 



y = const, fur kritische Druckverhaltnisse 



(2.2) 



bestimmt. 
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[0027] An der Drosselstelle, d.h. an der Drosselklappe auftretende Stromungsverluste werden Ober die geeignete 
Wahl von A^££) be rucks ichtigt. Aus stationaren Messungen kann bei bekannten Drucken vor und hinterder Drosselst- 
elle und bekanntem Massenstrom durch die Drosselstelle eine Zuordnung zwischen dem vom Drosselklappenstel- 
lungsfuhler 1 4 ermittelten Drosselklappenwinkel und dem entsprechendem reduzierten Querschnitt A/^gp angegeben 
werden. 

[0028] Wird der Luftmassenstrom an der Drosselklappe durch die Beziehung (2.2) beschrieben, so entsteht ein 
komplizierter Algorithmus zur numerisch richtigen Losung der Differentialgleichung (2.1 ). Zur Reduktion des Rechen- 
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aufwandes wird die DurchfluQfunktion \^ durch einen Polygonzug approximiert. 

[0029] Figur 2 zeigt den Verlauf der DurchfluBfunktion \|/ und das darauf angewandte Approximationsprinzip. Inner- 
halb eines Abschnitles i (I = 1...k) wird die DurchfluBfuntion \|/ dutch eine Gerade dargestelll. Mil einer vertretbaren 
Anzahl von Geradenabschnitten kann damit eine gute Approximation erreicht warden. Durch einen solchen Ansatz 
s kann die Gleichung (2.2) zur Berechnung des Massenstronnes an der Drosselklappe durch die Beziehung 
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(2.3) 



furl = (1...k) 
IS approximiert werden. 

[0030] In dieser Form beschreibt nrij die Steigung und n^ das Absolutglied des jeweiligen Geradenabschnlttes. Die 
Werte fur die Steigung und tOr das Absolutglied werden In Tabellen als Funktion des Verhaltnisses Saugrohrdruck zu 
Umgebungsdruck abgelegt. a 

[0031] Auf der Abfeisse von FIgur 2 1st dabei das Druckverhaltnis-^ und auf der Ordinate der Funktlonswert (0 - 
20 0.3) der DurchflulBfunktion v|/ aufgelragen. ^ 
[0032] Fur Druckverhaltnisse 
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ist v = konstant, d.h. daB der DurchfluB an der Drcsselstelle nur noch vom Querschnitt abhangig ist und nicht mehr 
von den Druckverhaltnissen. Die in die jeweiligen Zylinder der Brennkraftmaschine einstromende Luftmasse laBt sich 
analytisch nur schwer bestimmen. da sie stark vom Ladungswechsel abhangt. Die Fullung der Zylinder wird weitest- 
gehend durch den Saugrohrdruck. die Drehzahl und durch die Ventilsteuerzeiten bestimmt. 

[0033] Zur moglichst genauen Berechnung des Massenstroms in den jeweiligen Zylinder mz>, ist deshalb einerseits 
die Beschreibung der Verhaltnlsse im Ansaugtrakt der Brennkraftnnaschine mittels partieller Differentialgleichungen 
und andererserts die Berechnung des Massenstromes am EinlaBventil nach der DurchfluBgleichung als erforderliche 
Randbedingung notwendig. Erst dieser komplizierte Ansatz gestattet die Berucksichtigung dynamischer Nachladeef- 
fekte. die von der Drehzahl, der Saugrohrgeometrie, der Zylinderzahl sowie den Ventilsteuerzeiten maBgeblich beeln- 
fluBt werden. 

[0034] Da eine Berechnung nach oben genanntem Ansatz in der elektronlschen Steuerungseinrichtung der Brenn- 
kraftmaschine n^cht realisierbar ist, geht eine mogliche Naherung von einem einfachen Zusammenhang zwischen 
Saugrohrdruck Pq und Zyindermassenstrom Iz^ aus. Fur einen weiten Bereich der sinnvollen Ventilsteuerzeiten kann 
dafur in guter Naherung von einem linearen Ansatz der Form 



mzyl APPROX = yr-P5 + 7o (2.4) 
ausgegangen werden. 

[0035] Die Steigung 7^ und das Absolutglied 70 der Beziehung (2.4) sind dabei. unter Berucksichtigung aller wesenl- 
lichen EinfluBfaktoren Funktionen der Drehzahl, der Saugrohrgeometrie, der Zylinderzahl, der Ventilsteuerzeiten sowie 
der Temperatur der Luft im Saugrohr Tg. Die Abhangigkeit der Werte von 7^ und 7q von den EinfluBgroBen Drehzahl, 
Saugrohrgeometrie, Zylinderzahl und den Ventilsteuerzeiten und Ventilerhebungskurven kann dabei Ober stationare 
Messungen ermittelt werden. Ober diese Wertebestimmung wird ebenfails der EinluB von Schwingrohr- und/oder Re- 
sonanzsaugsystemen auf die von der Brennkraftmaschine angesaugte Luftmasse gut wiedergegeben. Die Werte von 
7^ und 7o sind in Kennfeldern der elektronlschen Motorsteuerungseinrichtung abgelegt. 

[0036] Als bestimmende GroBe zur Enmittlung der Motorlast wird der Saugrohrdruck Pg ausgewahlt. Mit Hilfe der 
Modell-Differentialgleichung soil diese GroBe moglichst exakt und schnell geschatzt werden. Die Schatzung von Pq 
erfordert die Losung der Gleichung (2.1). 

[0037] Mit den anhand der Formein (2.2) und (2.3) eingefOhrten Vereinfachungen kann (2.1) durch die Beziehung 
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(2.5) 



approximiert werden. Betrachtet man, entsprechend den Voraussetzungen zur Herleitung von Gleichung (2.1). die 
Tennperatur der Luft im Saugrohr Tg als eine langsam veranderliche Me6gr63e sowie A^^q als Eingangsgr63e, so 
Ia3t sich die nichtlineare Form der Differentialgleichung (2.1) durch die bilineare Gleichung (2.5) approximieren. 
[0038] Zur Losung der Gleichung (2.5) wird diese Beziehung In eine geeignete DIfferenzenglelchung ubergef uhrt. 
[0039] Als Kriterium zur Auswahl des geeigneten Dlfferenzenschemas konnen die folgenden prinzipiellen Anforde- 
rungen an die Losungseigenschaften der zur bildenden DIfferenzenglelchung fonmuliert werden: 



1 Das Differenzenschema muB auch unter extremen dynamlschen Anforderungen konsen^ativ sein, d.h. die Lo- 
sung der Differenzenglelchnung muB der Losung der Differentialgleichung entsprechen. 

2. die numerische Stabilltat muf3 zu Abtastzeiten, die den maximal moglichen Segmenlzeiten entsprechen, im 
gesamten Arbeitsbereich das Saugrohrdruckes garantiert sein. 



[0040] Forderung 1 ist durch einen impliziten Rechenalgorithmus erfullbar. Aufgrund der Approximation der nichtli- 
nearen Differentialgleichung (2.1) durch eine bilineare Gleichung ist das entstehende implizite Losungsschema ohne 
Einsatz iterativer Verfahren losbar, da die DIfferenzenglelchung in eine explizite Form uberfuhrt werden kann, 
[0041] Die zweite Forderung ist aufgrund der Konditionierung der Differentialgleichung (2. 1 ) und deren Approxima- 
tion (2.5) nur durch eine Rechenvorschrift zur Bildung der Differenzengleichung erfullbar, die absolut stabil arbeitet. 
Diese Verfahren werden auch als A-stabile Verfahren bezeichnet. Kennzeichnend fur diese A-Stabilitat ist die Eigen- 
schaft des Algorithmus. bei einem stabilen Ausgangsproblem fur beliebige Werte der Abtastzeit, d.h. Segmentzeit 
numerisch stabil zu sein. Eine mogliche Rechenvorschrift zur numerischen Losung von Differentialgleichungen. die 
beiden Forderungen gerecht wird, ist die Trapezregel. 

[0042] Die durch Anwendung der Trapezregel entstehende Differenzengleichung lautet im vorliegenden Fall 
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p;[;v-i]+p;[A^] 
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fur N (1. 



(2.6) 



definiert. 

[0043] Wird diese Vorschrift auf (2.5) angewandt, so ergibt sich die Beziehung 
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[N-\] + ^ Ps [N-]] 
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i 2k 1 
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2 y 



2k 1 



K-1 Rj_Ts 



1 Hl-Ts 



^RED 



I 2k 1 



fur N = und i = (l.-.k) (2.7) 

A 

zur Berechnung des Saugrohrdruckes PsW MaB fur die Motorlast. 

[0044] [N] bedeutet dabei das aktuelle Segment bzw. der aktuelle Rechenschritt. [A/+1 ] das nachstfolgende Segment 
bzw. def nachstfolgende Rechenschritt. 

[0045] Im folgenden wird die Berechnung des^aktuellen und pradizierten Lastsignales beschrieben. 
[0046] Aus dem berechneten Saugrohrdruck Pq kann der Luftmassenstrom der in die Zylinder einstromt, durch 
die Beziehung (2.4) enmittelt werden. Wendet man einen einfachen Integrationsalgorithmus an, so erhalt man fur die 
wahrend eines Ansaugtaktes von der Brennkraftmaschine angesaugte Luftmasse die Beziehung 



30 



mzv/[A']-^-^mzv/ [A^-l] + mzyt [N] 



3S 



fiir N - (1. . jQP) 



(2.8) 



[0047] Dabei wird davon ausgegangen. da3 der Anfangswert der Lastgr63e null ist. Fur die segmentsynchrone La- 
sterfassung sinkt mit steigender Drehzahl die Segmentzeit, wahrend die Segmentanzahl. urn die eine Kraftstoffvorla- 
gerung vorgenommen wird, steigen muB. Aus diesem Grund Ist es erforderlich, die Pradiktion des Lastsignals fur einen 
veranderlrchen Pradiktionshorizont H, d.h. fur eine bestimmte, in erster Linie drehzahlabhanglge Anzahl H von Seg- 
menten, auszulegen. BerOcksichtigt man diesen veranderlichen Pradiktionshorizont H, so kann Glerchung (2.8) In der 
Form 



so 



mzyi 
ftar N 



(1. . OO) 



(2,9) 



geschrieben werden. 

[0048] Fur die weiteren Uberlegungen wird davon ausgegangen. dafi sich die Segmentzeit Ty^ und die Parameter 
Yi und Yo de** Beziehung (2.4). die zur Bestimmung des Massenstromes i.^^ aus dem Saugrohrdruck Pq erforderlich 

sind, Ober die Pradiktionszeit nicht andem. Unter dieser^Voraussetzung wird die Pradiktion eines Wertes fur i.A*[A/+H] 
durch die Pradiktion des entsprechenden Druckwertes PslN+H\ erreicht. Dadurch nimmt die Gleichung (2.9) die Form 
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fur N = (1 . . ,<P) 



(2.10) 



an. 



10 [0049] Da bei dem beschriebenen Verfahren die zeitliche Anderung des Saugrohrdruckes in analytischer Form 
vorliegt, wird im folgenden die Pradiktion des Druckwertes [N+H\ durch H-fache Anwendung der Trapezregel er- 
reicht. In diesem Fall erhalt man die Bezrehung 
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Fs [N^H]^Ps [N] +^'H 



20 



fur N « (1- .90 ) 



(2.11) 
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[0050] Bestimmt man den Druck Pg[A/+H-1] in analoger Weise, so kann fur das pradizierte Lastsignal die Gleichung 



mzyl 



fur N - {1. .<X» 



(2-12) 
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angegeben werden. 

[0051] Wahit man fur den Pradiktionshorizont H Werte in der GroBen-ordnung von 1 .3 Segmenten, so kann mit der 
Fomriel (2.12) ein gut pradizlertes Lastsignal erhalten werden. 

[0052] Im folgenden wird das Prinzip des Modellabgleichs fur luftmassen- und saugrohrdruckgefuhrte Motorsteue- 
rungssysteme erklart. 

[0053] Bedingt durch den Einsatz von Motoren mit variabler Ventilsteuerung und/oder veranderlicher Saugrohrgeo- 
metrie, durch Fertigungstoleranzen und Alterungserscheinungen, sowie durch Temperatureinflusse sind die Werte von 
7i und Yo nnit einer gewissen Unsicherheit behaftet. Die Parameter der Gleichung zur Bestimmung des Massenstromes 
in den Zyllndem sind, wie oben beschrieben, Funktionen vielfaltiger EinfluBgr63en, von denen nur die wichtigsten 
erfafBt werden konnen. 

[0054] Bei der Berechnung des Massenstromes an der Drosselklappe wirken sich Mef3fehler bei der Erfassung des 
Drosselklappenwinkels und Approximationsfehler bei der Polygonzugapproximation der DurchfluBfunktlon \f auf die 
ModellgroGen aus. Besonders bei kleinen Drosselklappenwinkein ist die Syslemempfindlichkeit gegenuber erstge- 
nannten Fehlern besonders hoch. Daraus ergibl sich, daB kleine Anderungen der Drosselklappenstellung einen gra- 
vierenden EinluB auf Massenstrom bzw. Saugrohrdruck haben. Um die Wirkung dieser Einflusse zu reduzieren, wird 
im folgenden ein Verfahren vorgeschlagen, das es gestattet. bestimmte GroBen. die EinfluB auf die Modellrechnung 
haben, so zu korrigieren, daB eine genauigkeitsverbessernde Modellanpassung fur stationaren und instationaren Mo- 
torbetrieb durchgef uhrt werden kann. 

[0055] Die Anpassung wesentlicher Parameter des Modells zur Bestimmung der LastgroBe der Brennkraftmaschine 
erfolgt durch die Korreklur des aus dem gemessenen Drosselklappenwinkel bestimmten reduzierten Querschnitts A^^q 
durch die KorrekturgroBe AA^^ 

[0056] Die EingangsgroBe zur korrlgierten Saugrohrdruckberechnung kf^^p wird damit durch die Beziehung 
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^REDKORR - ^RED ^^RED (3-11) 

beschrieben. 

[0057] In der Gleichung (2.2) und nachfolgenden Formein wird dann A/^^^ durch f^^EDKORR ersetzt. Zur Verbesse- 
rung des Folgeverhaltens des Regelkreises wird der aus dem MeBwert des Drosselklappenwinkels abgeleitete redu- 
zierte Drosselklappenquerschnitt A^fijin die Modellrechnung einbezogen. Die KorrekturgroBe AA^^gu wird durch Rea- 
iisierung eines Modellregelkreises gebildet. 

[0058] Fur luftmassengef uhrte Motorsteuerungssysteme ist der mittels des Luftmassennnessers an der Drosselklap- 
pe gemessene Luftmassenstrom nio/c_LA^iw cJ'© FuhrungsgroBe dieses Regelkreises, wahrend fur saugrohrdruckge- 
f uhrte Systeme der gemessene Saugrohrdruck P^als FuhrungsgroBe genutzt wird. Uber eine Folgeregelung wird der 
Wertvon AAflPoSobestimmt, daBdie Regelabweichung zwischen FuhrungsgroBe undderensprechenden RegeigroBe 
minimiert wird. 

[0059] Um auch im dynamischen Betrieb Genauigkeitsverbesserungen mit der genannten Methode zu erreichen, 
muB die MeBwerterfassung der FuhrungsgroBe moglichst exakt nachgebildet werden. In den nneisten Fallen sind dabei 
das dynamische Verhalten des Sensors, d.h. entweder des Luftmassenmessers Oder des Saugrohrdrucksensors und 
eine nachfolgend durchgefuhrte Mittelwertbildung zu berucksichtigen. 

[0060] Das dynamische Verhalten des jeweiligen Sensors kann in erster Naherung als ein System erster Ordnung 
mit eventuell arbeitspunktabhangigen Verzogerungszeiten modelliert werden. Im Falle eines luftmassengef Oh rten 
Systems lautet eine mogtiche Gleichung zur Beschreibung des Sensorverhaltens 



2B 



mDK_iJ4M [N]^e^' ' inoK [A^-l]+ VinDK^iMMlN ^\] (3.12) 



[0061] Eine GroBe. die beim gewahlt^n Ansatz einen wesentlichen EinfluBauf den maximal moglichen Massenstrom 
inz^ besitzt, ist der Umgebungsdruck Pfj. Aus diesem Grund kann nicht von einem konslanten Wert dieser GroBe 
ausgegangen werden, sondern es erfolgt eine Anpassung in der nachfolgend beschriebenen Art und Weise. 
[0062] Der Wert des Umgebungsdruckes wird verandert, ^enn der Betrag der KorrekturgroBe AA^^^ 
stimmte Schwelle uberschreitet oder wenn das Druckverhaltnis-pS groBer als eine wahlbare Konstante ist. Damit wird 
gewahrleistet, daB sowohl im Teil- als auch im Vollastbereich emU Umgebungsdruckanpassung erfolgen kann. 
[0063] Im folgenden wird ein Model labgleich fur luftmassengefuhrte Motorsteuerungssysteme erklart. Fur dieses 
System kann die in Figur 3 dargestellte Modellstruktur angegeben werden. 

[0064] Der Drosselklappenstellungsfuhler 14 (Figur 1) liefert ein dem Offnungsgrad der Drosselklappe 11 entspre- 
chendes Signal, z.B. einen Drosselklappenoffnungswinkel. In einem Kennfeld der elektronischen Motorsteuerungs- 
einrk^htung sind zu verschiedenen Werten dieses Drosselklappenoffnungswinkels zugehorige Werte fur den reduzler- 
ten Querschnitt der Drosselklappe kf^ED abgespeichert. Diese Zuordnung wird durch den Block "statisches Modell" in 
Figur 3 und In Figur 4 reprasentiert. Das Teilsystem "Saugrohrmodell" in den Figuren 3 und 4 reprasentiert das durch 
(2.7) beschriebene Verhalten. FuhrungsgroBe dieses Modellregelkreises ist der MeBwert des uber ein Segment ge- 
mittelten Luftmassenstromes an der Drosselklappe Wird als Regler in diesem Modellregelkreis ein Pl-Regler 

eingesetzt. so ist die bleibende Regelabweichung null, d.h. ModellgroBe und MeBgroBe des Luftmassenstromes an 
der Drosselklappe sind identisch. Die Pulsationserscheinungen des Luftmassenstromes an der Drosselklappe. die vor 
allem bei 4-Zylindermotoren zu beobachten sind, fuhren bei betragsbildenden Luftmassenmessern zu erheblichen 
positiven MeBfehlern und somit zu einer stark fehlerbehafteten FuhrungsgroBe. Durch eine Abschaltung des Reglers, 
d.h. einer Verkleinerung der Reglerparameter kann zum gesteuerten modellgestOtzten Betrieb ubergegangen werden. 
Bereiche, in denen die genannten Pulsatk>nen auftreten. konnen somit mit dem selben Verfahren unter Berucksichti- 
gung dynamischer Zusammenhange behandelt werden, wie diejenlgen Bereiche, In denen eine nahezu ungestorte 
FuhrungsgroBe vorliegt. Im Gegensatz zu Verfahren, die relevante MeBwerte nur in statlonaren Betriebspunkten be- 
rucksichtigen, bleibt das beschriebene System nahezu uneingeschrankt arbeitsfahig. Bei Ausfall des Luftmassensi- 
gnals oder des Signals des Drosselklappensteltungsf uhlers ist das vorgestellte System in der Lage, ein entsprechendes 
Ersatzsignal zu bilden. Bei Ausfall der FuhrungsgroBe muB der gesteuerte Betrieb realisiert werden, wahrend im an- 
deren Fall der geregelte Betrieb die kaum beeintrachtigte Funktionsfahigkeit des Systems garantiert. 
[0065] Der Block "Saugrohrmodell" reprasentiert die Verhaltnisse wie sie anhand der Gleichung (2.7) beschrieben 
sind und hat demzufolge als AusgangsgroBe die ModellgroBe P^sowie die zeitliche Ableitung >^ und die GroBe ^m-. 
Nach der Modellierung des Sensorubertragungsverhaltens d.h. des Ubertragungsverhattens des Luftmassenmessers 
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und der Abtastung wird die Modellgr63e einer Mittelung unterzogen, so daB die gemittelte GroBe und der 
vom Luftmassenmesser gemessene durchschnittliche Luftmassenstrom einem Vergleicher zugefuhrt werden 
konnen. Die Differenz beider Signale bewirkt eine Anderung AA^f^pdes reduzierlen Strom ungsquer-schnlttes Afj^ 
so da3 stationar und instationar ein Mo-dellabgleich erfolgen kann. 

[0066] Fur saugrohrdruckgefuhrte Motorsteuerungssysteme wird die in Figur 4 dargestellte Modellstruktur angege- 
ben, wobei gleiche Blocke wie in Figur 3 gleiche Bezeichnungen tragen. Ebenso wie bei dem luftmassengetuhrten 
Motorsteuerungssystem reprasentiert das Teilsystem "Saugrohrnnodell", das durch die Differenzengleichung (2.7) be- 
schriebene Verhalten. FuhrungsgroBe dieses Modellregelkreises 1st der MeBwert des uber ein Segment gemittelten 
Saugrohrdruckes Pg_^ Wird ebenfalls wie in Figur 3 ein Pl-Regler eingesetzt, so 1st Im stationaren Fall der MeBwert 
des Druckes im Saugrohr mit der ModellgroBor, . identisch. Wie oben beschrieben, bleibt auch das vorliegende 
System nahezu uneingeschrankt arbeitsfahig, da bei Ausfall des Saugrohrdrucksignales oder des MeBwertes fur den 
Drosselklappenwinkel ein entsprechendes Ersatzsignal gebildet werden kann. 

[0067] Die durch das Saugrohrmodell erhaltenen ModellgroBen P3, h werden einem Block "Pradiktion" zugefuhrt. 
Da mit den Modellen auch die Druckanderungen im Saugrohr berechnet werden, konnen diese Druckanderungen dazu 
verwendet werden. den zukunftigen Druckverlauf im Saugrohr und damit die Zylinderluftmasse fur das nachste [A/+1] 
Oder f urdie nachsten Segmente [A/+H] zu schatzen. Die GroBe /772y/bzw. die GroBe /f?^y/[A/+1] dienen dann zur exakten 
Berechnung der Einspritzzeit, wahrend derer Kraflstoff eingespritzt wird. 



Patentanspruche 

1. Vedahren zum Bestimmen der einstromenden Luftmasse in den bzw. die Zylinder einer Brennkraftmaschine mit 

einem Ansaugsystem, das ein Saugrohr (10) und eine darin angeordnete Drosselklappe (11), sowie einen 
den Offnungsgrad der Drosselklappe (19) erfassendcn Drosselklappensteilungsfuhler (14) aufweist, 
einem, ein Lastsignal ( ; Ps_s) Brennkraftmaschine erzeugenden Sensor (12; 13) 
einer elektrischen Steuerungseinrichtung, die auf der Grundlage des gemessenen Lastsignals ( ; q) 
und der Drehzahl der Brennkraftmaschine eine Grundeinspritzzeit berechnet, wobel 

die Verhaltnisse im Ansaugsystem mittels eines Saugrohrfullungsmodells nachgebildet werden, wobei als 
EingangsgroBen des Modells der Offnungsgrad der Drosselklappe (11), der Umgebungsdruck {PJ) und die 
Ventilstellung reprasentierende Parameter herangezogen werden, 
- eine ModellgroBe fur den Luftmassenstrom (ir,*. ) an der Drosselklappe (11 ) mit Hitfe der DurchfluBgleichung 
Idealer Gase durch Drosselstellen beschrieben wird, 

eine ModellgroBe^fur den Luftmassenstrom (mz.) in den bzw. in die Zylinder (17) als lineare Funktion des 

Saugrohrdrucks (P5) durch eine Massenbilanz der Luftmassenstrome (i,^ , ^2.*) beschrieben wird 

diese ModellgroBen uber eine Differentialgleichung verknupft werden, daraus als bestimmende GroBe zur 

Ermittlung der tatsachlichen Last der Brennkraftmaschine der Saugrohrdruck (P^) berechnet wird und 

aus dem linearen Zusammenhang zwischen berechnetem Saugrohrdruck (P5) und der ModellgroBe fur den 

Luftmassenstrom {^z^) in den bzw. in die Zylinder (17) durch Integratbn die in den bzw. in die Zylinder (17) 

einstromende Luftmasse (m^iy,) erhalten wird. 

2. Verfahren nach Anspruch 1 , 
dadurch gekennzeichnet, 

daB das vom Lastsensor (12; 13) gemessene Lastsignal ( l,o*._ua, ; P5f_s)zur Korrektur und damit zum Abgleich der 
ModellgroBen (i.?^) in einem geschtossenen Rgelkreis herangezogen wird, wobei das Lastsignal ( ; Ps_s) 
als FQhrungsgroBe des Regelkrelses dient. 

3. Verfahren nach Anspruch 2, 
dadurch gekennzeichnet, 

daB der Abgleich im stationaren und/oder instationaren Betrleb der Brennkraftmaschine durchgefuhrt wird und 
dabei das Ubertragungsverhalten des l^stsensors (12; 13) berucksichtigt wird. 

4. Verfahren nach Anspruch 2, 
dadurch gekennzeichnet, 

daB jedem gemessenen Wert des Drosselklappenoffnungsgrades ein Wert eines reuzierten Querschnitts der Dros- 
selklappe (Af^EQ) zugeordnet ist und der Abgleich der ModellgroBen durch Korrektur des reduzierten Querschnitts 
i^REo) durch eine KorrekturgroBe {^r^d) derart erfolgt, daB die Regelabwek:hung zwischen FQhrungsgroBe und 
entsprechender ModellgroBe minimiert wird. 
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5. Verfahren nach Anspruch 4, 
dadurch gekennzeichnet, 

daB der reduzierte Querschniti {Af^^p) aus slalionaren Messungen am Motorprufstand ermittelt wird und in einem 
Kennfeld eines Speichers der elektrischen Steuerungseinrichtung abgelegt ist. 

5 

6. Verfahren nach Anspruch 1 , 
dadurch gekennzeichnet, 

da3 bei der Darstellung der ModellgroBe fur den Luftmassenstrom (^o*) an der Drosselklappe (11) eine in der 
DurchfluBgletchung (Gleichung 2.2) vorhandene DurchflufBfunktion (y) in einzelne Abschnltte (i = 1...k) unterteilt 
10 wird und diese Abschnltte durch Geradenabschnitte angenahert werden, wobei fur die Steigung (nr\j) und fur das 

Absolutglied {n^) der jeweiligen Geradenabschnitte als Funktion des Verhallnisses von Saugrohrdruck (P^) und 
Umgebungsdruck (PJj bestimmt sind und in einem Kennfeld abgelegt sind. 

7. Verfahren nach Anspruch 1 . 
IS dadurch gekennzelchnet, 

daB die Steigung (y^ ) und das Absolutglied (yq) der linearen Funktion fur die ModellgroBe fur den Luftnnassenstrom 
in den bzw. ind die Zylinder ( my.. ) abhangig von nnindestens einem der Parameter, Drehzahl der Brennkrattmaschi- 
ne, Zylinderzahl, Saugrohrgeometrie, Temperatur der Luft (Tq) Im Ansaugrohr (10) und Ventllsteuerzeichen fest- 
gelegt sind. 
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8. Verfahren nach Anspruch 7, 
dadurch gekennzeichnet, 

daB die Parameter duch stationare Messungen am Motorprufstand ermittelt werden und in Kennfelder abgelegt 
sind. 

9. Verfahren nach Anspruch 1. 
dadurch gekennzelchnet, 

daB die in den Zylinder einstromende Luftmasse (m^y/) durch die Beziehung 



berechnet wird, mit 

T^: Abtastzeit oder Segmentzeit 
40 [A/] ModellgroBe des Luftmassenstromes wahrend des aktuellen Abtastschrittes Oder Segments 

ii»^ [A(-1]ModellgroBe des Luftmassenstromes wahrend des vergangenen Abtastschrittes oder Segmentes. 

10. Verfahren nach Anspruch 1, 
dadurch gekennzelchnet, 

daB die in den bzw. in die Zylinder einstromende Luttmasse (rhz^ fur einen, bezuglrch der aktuellen Lasterfassung 
zum Abtastzeitpunkt [N\ in der Zukuntt liegenden, bestimmten Pradiktionshorizont (H) geschatzt wird durch Schat- 
zung des entsprechenden Druckwertes nach folgender Beziehung: 



so 



ss mit 



T^: Abtastzeit oder Segmentzeit 

H: Pradiktionshorizont. Anzahl der In der Zukunft liegenden Abtastschrltte 
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7^ : Steigung der linearen Gleichung 
7o- Absolutglied zur Bestimmung von 
N: aktuelier Abtastschritt 

11. Verfahren nach Anspruch 10, 
dadurch gekennzeichnet, 

daB die Anzahl (H) von Segmenten, fur die das Lastsignal fur die Zukunft geschatzt werden soil, drehzahlabhangig 
f estgelegt ist. 

Claims 

1 . Metliod for determining the air mass flowing into the cylinder or cylinders of an internal combustion engine with 

^5 - an induction system which comprises an Induction manifold (1 0) and a throttle valve (11 ) arranged within it as 

well as a throttle valve position sensor (14) which detects the degree of opening of the throttle valve (19), 

a sensor (12; 13) which produces a load signal ( mtiK iUM s) the internal combustion engine, 

20 - an electrical control device, which calculates a basic injection time on the basis of the measured load signal 

( iw.tw ; Ps s) and the speed of the internal combustion engine, whereby 

the conditions in the induction system are simulated by an Induction manifold charging model, wherein pa- 
rameters representing the degree of opening of the throttle valve (11 ), the ambient pressure (P^) and the valve 
25 position are used as input variables of the model, 

a model variable for the air mass flow (^^^ ) at the throttle valve (11 ) is described by means of the flow-through 
equation of ideal gases through throttling points, 

30 - a model variable for the air mass flow (Lw ) into the cylinder or cylinders (1 7) is described as a linear function 

of the induction nnanifold pressure (Ps) by a mass balance of the air mass flows m«w). 

A 

these model variables are linked by a differential equation, the inductbn manifold pressure {Ps) being calcu- 
lated therefrom as a defining variable for determining the actual load of the internal combustion engine, and 

35 

the air mass (wzw) flowing into the cylinder or cylinders (17) is obtained from the linear relationship between 
calculated induction manifold pressure (Pq) and the model variable for the air mass flow (m^) in the cylinder 
or cylinders (17) by integration. 

40 2. Method in accordance with claim 1, characterized in that the load signal (^w.iw*i; Ps s) Pleasured by the load 
sensor (12; 13) is used to correct and therefore to adjust the model variables {miy,) in a closed control circuit, the 
load signal ( i,M_tM»; T^s_s) serving as reference variable of the control circuit. 

3. Method in accordance with claim 2, characterized in that the adjustment is made in the steady-state and/or transient 
45 mode of the internal combustion engine, allowing for the transmission behaviour of the load sensor (12; 1 3). 

4. Method in accordance with claim 2, characterized in that a value of a reduced cross-section of the throttle valve 
(^red) allocated to each measured variable of the degree of opening of the throttle valve and the model variables 
are adjusted by correcting the reduced cross-section (^red) ^ correction variable (A^f^^o) such that the deviation 

so between reference variable and corresponding model variable is minimized. 

5. Method in accordance with claim 4. characterized in that the reduced cross-section (A^^^q) is determined from 
steady-state measurements on the engine test stand and is entered into a map in a memory in the electric control 
device. 



ss 



Method in accordance with claim 1 , characterized in that, when portraying the model variable for the air mass flow 
(moK) at the throttle valve (11), a flow-through function (\|/) present in the flow-through equation (equation 2.2) is 
divided into individual sections (i = 1...k) and these sections are approximated by straight line sections, the re- 
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spective straight line sections being detornnined as a function of the ratio of induction manifold pressure (P^) to 
ambient pressure (P^) for the gradient (mj) and for the absolute term (nj) and being entered into a map. 

7. Method in accordance with claim 1 , characterized in that the gradient (yi) and the absolute term (yq) of the linear 
s function for the model variable for the air mass flow into the cylinder or cylinders ( «w ) are established as a function 

of at least one of the parameters of speed of the internal combustion engine, number of cylinders, induction manifold 
geometry, temperature of the air (Tg) in the induction manifold (10) and valve control symbol, 

8. Method in accordance with claim 7, characterized in that the parameters are determined by steady-state meas- 
10 urements on the engine test stand and are entered into maps. 

9. Method in accordance with claim 1 , characterized in that the air mass (m^jj flowing into the cylinders is calculated 
by the equation 



IS 



20 



25 



30 



where 

is the sampling interval or segment time 

m2yi[N\ is the model variable of the air nnass flow during the current sampling step or segment 

mz»t [AM] is the model variable of the air mass flow during the previous sampling step or segment. 

10. Method in accordance with claim 1 , characterized in that the air mass (m^y/) flowing into the cylinder or cylinders 
is estimated for a defined prediction horizon (H) located in the future with respect to the current load detection at 
the moment of sampling [Afl by estimating the corresponding pressure value in accordance with the following 
equation: 
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where 

Ta: 
H: 

Yi : 

7o: 
N: 



sampling interval or segment time 

prediction horizon, number of sampling steps located in the future 
gradient of the linear equation 
absolute term for determining m^.. 
current sampling step 



11. Method in accordance with claim 10. characterized in that the number (H) of segments for which the load signal 
is to be estimated for the future is established as a function of the speed. 



ss Revendications 

1. ProcddS de determination du d6bit massique d'air introduit dans le ou les cylindres d'un moteur d combustion 
interne, comprenant : 
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un systeme d'admission comportant un cotlecteur d'admission (10) et un papillon des gaz (11) dispose dans 
celui-ci, ainsi qu'un dispositif capteur de position de papillon des gaz (14) detectant le degr6 d'ouverture du 
papillon des gaz (19), 

un capteur (12 ; 13) produisant un signal de charge ( ; Ps_s) moteur a combustion interne, 
un dispositif de commando electronique qui calcule une duree d'injection de base k partir du signal de charge 
mesure ( ^r...»« ; Ps s) et de la Vitesse de rotation du moteur a combustion interne, 

selon lequel : 

les conditions dans le systeme d'admission sont simul6es au moyen d'un modete de remplissage de coilecteur 
d'admission, le degre d'ouverture du papillon des gaz (11), la pression amblante (PJ) et des parametres re- 
presentant la position de soupapes etant utilises en tant que grandeurs d'entree du modele, 
une grandeur de modele correspondant au d6bit massique d'air (^o* ) a I'endroit du papillon des gaz (11) est 
decrite au moyen de Inequation d'ecoulement des gaz parfaits a travers des emplacements d'etranglement, 
une grandeur de modele correspondant au debit massique d'aIr (mz^) dans le ou les cylindres (17) est decrite 
comme fonction lineaire de la pression de coilecteur d'admission (P^) au moyen d'un biian des debits massi- 
ques d'air (moK , "z.i ), 

ces grandeurs de modele sont iiees au moyen d'une Equation differentlelle et, a partir de celle-ci, la pression 
de coilecteur d'admission (P^) est calcul6e comme grandeur d6cisive pour la determination de la charge ef- 
fective du moteur a combustion inlerne el ^ 
k partir de la relation lineaire se presentant entre la pression de coilecteur d'admission calculee (P3) et de la 
grandeur de modele correspondant au debit massique d'air («z« ) dans le ou les cylindres (17), le debit d'air 
(^Zy/) introduit dans le ou les cylindres (17) est obtenu par integration. 

Procede suivant la revendication 1 . caractcrisc on ce que le signal de charge ( 'mo. j.uu ; Ps_s) f^iesure par le capteur 
de charge (12 ; 13) est utilise dans un circuit de regulation terme pour la correction et done pour la compensation 
des grandeurs de modele (^... ) le signal de charge ( ; P^ s) sen^ant de grandeur de pilotage du circuit de 
regulation. 

Procede suivant la revendication 2. caracterise en ce que la compensation est effectu6e en regime permanent et/ 
ou en regime non permanent du nnoteur a combustion interne et qu'il est de ce fait tenu compte du comportement 
de transmission du capteur de charge (12:1 3). 

Procede suivant la revendication 2, caracterise en ce qu'une valeur d'une section transversale r^duite (Afj^Q) du 
papillon des gaz est associee a chaque valeur mesuree du degre d'ouverture de papillon des gaz et que la com- 
pensation des grandeurs de modele par correction de la section transversale reduite (A^ed) moyen d'une 
grandeur de correction (^red) ®st effectuee de fagon telle que I'ecart de regulation entre la grandeur de pilotage 
et la grandeur de module correspondante est rendu minimal, 

Procede suivant la revendication 4, caracterise en ce que la section transversale reduite (Ared) ®st determinee 
a partir de mesures effectuees en regime permanent sur banc d'essais de moteur et est rangee en memoire dans 
une table caracteristique d'une memoire du dispositif de commando eiectrique. 

Procede suivant la revendication 1, caracterise en ce que, pour la representation de la grandeur de modele cor- 
respondant au debit massique d'air (i„oK ) a I'endroit du papillon des gaz (11 ), une tonction d'ecoulement (\^) existant 
dans {'equation d'ecoulement (equation 2.2) est divisee en segments individuels (i = 1 . . . k) et ces segments font 
I'objet d'une approximation au moyen de segments de droite, la pente (rc\) et I'abscisse k I'origine (n^) des segments 
de droite respectifs etant determines en fonction du rapport de la pression de coilecteur d'admission (P^ el de la 
pression ambiante (Py) et etant rangees en menrraire dans une table caracteristique. 

Procede suivant la revendication 1. caracterise en ce que la pente (7^) et I'abscisse a I'origine (yo) de la fonction 
lineaire correspondant a la grandeur de modele prevue pour le debit massique d'air (m^) dans te ou les cylindres 
sont determinees en fonction d'au moins I'un des parametres que sont la Vitesse de rotation du moteur ^ combustion 
interne, le nombre de cylindres, les parametres geometriques du coilecteur d'admission, ia temperature (T^) de 
I'air dans le coilecteur d'admission (10) et les durees de commande de soupapes. 

Procede suivant la revendication 7, caracterise en ce que les parametres sont determines au moyen de mesures 
effectuees en regime permanent sur banc d'essais de moteur et sont ranges en memoire dans des tables carac- 
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tdristiques. 

9. Procede suivant la revendication 1. caracterise en ce que te debit massique d'air {mzyd penetrant dans le cylindre 
est calculi au moyen de la relation. 



mzyi[N]^^^\ m2yi[N-\]^mzy{[N] 



ddsignant la dur6e d'analyse par echantillonnage ou segment de temps, 
m^.ilN] la grandeur de modele du debit massique d'air pendant le pas d'analyse ou segment actuel. 

^zw (A/-1 ] la grandeur de modele du debit massique d'air pendant le pas d'analyse ou segment pass6, 

IS 10, Procede suivant la revendication 1. caracterise en ce que le d^bit massique d'air (m^y/) penetrant dans (e ou les 
cylindres est estim6 pour un horizon de prediction (H) determine, situe dans le futur par rapport a la determination 
de charge actuelle ^ I'instant d'analyse par echantillonnage [Afl, par estimation de la valeur de pression corres- 
pondante conformement ^ la relatbn suivante : 



avec 



: duree d'analyse par echantillonnage ou segment de temps, 



H 

Yi 

N 



horizon de prediction, nombre des pas d'analyse par Echantillonnage se trouvant dans le futur. 

pente de I'Equation lin^aire, 

abscisse h I'orlgine servant a la detemnination de 

pas d'analyse par echantillonnage actuel. 



11. PtocM6 suivant la revendication 10, caract6ris6 en ce que le nombre (H) de segments pour lesquels le signal de 
3S charge doit etre estlm6 dans le futur est determine en fonction de la Vitesse de rotation. 
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